Draft version December 12, 2012 

Preprint typeset using I^'T^]X style emulatcapj v. 8/17/12 



THE STRUCTURE OF THE BROAD LINE REGION IN AGN: 
1. RECONSTRUCTED VELOCITY-DELAY MAPS 

C. J. Grier\ B. M. Peterson^ Keith Horne'*, M. C. Bentz'', R. W. Pogge^^^, K. D. Denney^, G. De Rosa\ 
Paul Martini^ ^ C. S. Kochanek^-^ Y. Zu\ B. Shappee\ R. Siverd'', T. G. Beatty\ S. G. Sergeev^ S. Kaspi**-'-', 
C. Araya Salvo\ J. C. Bird\ D. J. Bord"\ G. A. Borman^-", X. CHE^^ C. CHEN^^ S. A. Cohen", M. Dietrich\ 
V. T. Doroshenko'"'^^, Yu. S. Efimov^*, N. Free", I. Ginsburg", C. B. Henderson^ A. L. KING^^ K. Mogren^ 
M. MoLiNA^'^", A. M. Mosquera\ S. V. Nazarov'-", D. N. Okhmat^", O. Pejcha\ S. Rafter'-', J. C. Shields^'', 
J. Skowron^, D. M. Szczygiel\ M. Valluri", and J. L. van Saders^ 

Draft version December 12, 2012 

ABSTRACT 

We present velocity-resolved reverberation results for five active galactic nuclei. We recovered 
velocity-delay maps using the maximum-entropy method for four objects: Mrk 335, Mrk 1501, 3C 120, 
and PC 2130-)-099. For the fifth, Mrk 6, we were only able to measure mean time delays in different 
velocity bins of the H/3 emission line. The four velocity- delay maps show unique dynamical signatures 
for each object. For 3C 120, the Balmer lines show kinematic signatures consistent with both an 
inclined disk and infalling gas, but the He ii A4686 emission line is suggestive only of inflow. The 
Balmer lines in Mrk 335, Mrk 1501, and PC 2130-1-099 show signs of infalling gas, but the Heii 
emission in Mrk 335 is consistent with an inclined disk. We also see tentative evidence of combined 
virial motion and infalling gas from the velocity-binned analysis of Mrk 6. The maps for 3C 120 and 
Mrk 335 are two of the most clearly defined velocity-delay maps to date. These maps constitute a 
large increase in the number of objects for which we have resolved velocity-delay maps and provide 
evidence supporting the reliability of reverberation-based black hole mass measurements. 
Subject headings: galaxies: active — galaxies: nuclei — galaxies: Seyfert 



1. INTRODUCTION 

One of the continuing mysteries in the studies of active 
galactic nuclei (AGN) is the structure and kinematics of 
the broad line region (BLR). It is generally accepted that 
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AGN are powered by supermassive black holes, and the 
broad emission lines seen in Type 1 AGN are the re- 
sult of the photoionization of gas in the BLR. Under the 
assumption that this gas is in virial motion around the 
black hole, we can use it to measure the mass of the black 
hole itself. However, it is impossible to study the struc- 
ture of the BLR directly because the BLR is on the order 
of light-days in radius, which renders it spatially unre- 
solvable even in the nearest of galaxies with the largest 
diffraction-limited telescopes. The unknown BLR struc- 
ture introduces uncertainties in any black hole mass mea- 
surements derived from BLR gas kinematics. To learn 
about the BLR structure, we must rely on either rever - 
beration mapping techniques fe.g.. iPeterson et aI]|2004D 
or microlensing of gravitationally lensed quasars (e.g., 
iGuerras et al.l 120121) . Reverberation techniques use the 
time variability observed in the AGN continuum emis- 
sion and the sub sequent response of the gas in the BLR 
(jBlandford fc McKec 1982 : Peterson. 1993) . By monitor- 
ing the AGN spectra over a period of time, we can de- 
termine the distance of the emitting gas from the central 
source by measuring the time delay between variations 
in the continuum and the response of the emission lines. 
This time delay is due to the light-travel time between 
the continuum source and the BLR. 

The variations in the BLR emission-line flux, AL(y, i), 
are a convolution of the continuum flux varia- 
tions, AC(t), with the "transfer function", \l/(V,r) 
(jBlandford fc McKed [1981 ) . The transfer function de- 
pends on the temporal lag r between the line and con- 
tinuum emission, and the line-of-sight velocity V of the 
BLR gas. The relationship is expressed mathematically 
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AL{V,t) 



'^{V,T)AC{t -T)dT. 



(1) 



A main goal of reverberation mapping is to recover 
^'(y, r), called the "velocity-delay map" , which describes 
how the continuum flux variations give rise to BLR flux 
variations, and therefore contains information about the 
BLR geometry and kinematics. For example, a Keple- 
rian disk produces a velocity-symmetric structure with a 
wider/narrower range of velocities at smaller/larger de- 
lays. In contrast, radial flows give rise to asymmetric 
velocity structure, the signature of infall/outflow being 
smaller delays on the red/blue side of the velocity pro- 
file. Most previous reverberation studies have been lim- 
ited to measuring the mean time delay (t) for various 
emission lines. This allows us to estimate the mean ra- 
dius of the BLR (e.g..lPeterson et al.Jl2004l : IPennev et al.l 
l2010HBentz et al.ll2010aHGrier et al.ll20Ta) . but it reveals 
very little information about the detailed structure of the 
BLR. 

It has long been recognized that measuring emission- 
line time lags as a function of line-of-sight velocity 
provides a wa y to determine the gross kinematics of 
the B LR fe.g. . iBahcall et al. | [T97l iBlandford fc McKed 
1982t[B ochkar ev fc Antokhin 119821: iCapriotti et al.lll982l: 



Antokhin fc Bochkarevl[l983D . although it is only rela- 



tively recently that suitable data have become available 
to achieve this in practice. There were a number of early 
attempts to extract velocity-dependent time lags in emis- 
sion lines in order to probe the BLR kinematics, though 
these were generally frustrated by low data quality, time- 
sampling issues, or both. The first attempts that we 
ar e aware of to search for velocity-dependent lags were 
bvlGaskelll ([1988 1 for NGC 4151 and Koratkar & Gaskell 
(|1989| 1 for Fairall 9, in both cases for the C iv A1549 emis- 
sion line based on data obtained with the International 
Ultraviolet Explorer (a mean lag for F airall 9 had been 
obtained earlier bv iClavel et al.l [l989D . In both cases, 
detection of infall was claimed. We have, however, rean- 
alyzed both of these data sets using more modern meth- 
ods and find that the light curves are far too poorly 
sampl ed to support eithe r claim of a detection of infall. 
Later. lClavel et al.l 1)19900 revisited the case of NGC 4151 
with improved time sampling. To within the accuracy of 
their measurement, they found the red and blue wings of 
C IV to vary simultaneously. However, their mean sam- 
pling interval of 3.4 days was only slightly smaller than 
the C IV mean response time, which would have made 
it difficult to detect the subtle velocity-dependent lags 
that have b een detected in optical lines in recent years ( 
IPennev et al. 2009; Bcntz et al. 2010a). 

The first well-sampled AGN emission-line light curves 
became available as a result of intensive monitoring of 
NGC 5548 with l UE (IClavel et all Il99l and ground- 
based telescopes ([Peterson et al.l Il991t [Dietrich et al.l 
1199^ in 1988-89. These data led to the first suc- 
cessful attempts to recov er one- dimensional delay maps, 
5'(t) (Home et al. 199l1 : iKrolik et al.l[T99l[ ). Other at- 
tempts we re also made to search for velocity-dependent 
lags (e.g., IClavell fT99ll for NGC 5548), but they gen- 
erally concluded that any veloc ity-dependent lags re- 
mained unresolved. However, iCrenshaw fc Blackwelll 
([T99Q1) argued that CiV profile variations during part 



of the 1988-89 campaign on NGC 5548 indicated in- 
fall of BLR gas, although an infall signature was 
not detected using any other subset of the cam- 
paign data and the lUE spectra were generally rather 
marginal for AGN emission-line profile studies. Other 
sources were searched for velocity-dependent emission- 
line resp onses, but these searches were also unsuccess- 
ful (e.g., iKollatschnv fc Dietrich! I1997D . Some limited 
succes s was achieved with both space-based obs erva- 
tions (lUlrich fc Hom e' '1996: 'Wanders et al.l 119951 ) and 
ground-based data (KoUatschnv 2003), but e xcept in 
the case of NGC 4151 ([Ulrich fc Horn j [19961 ) the re- 
sulting velocity-delay maps showed little clear struc- 
ture. There was also a growing understanding that 
the dramatic emission-line profile variations observed 
in some sources were not a reverberation e ffect (e.g., 
iWanders fc Peterson 1996: Scrge ev et al.ll2007| ). and that 
reverberation signals were generally going to be weak. 
Obtaining a high-fidelity velocity-delay map was go- 
ing to require bot h high-quality dat a and excellent 
temporal sampling ([Home et al.l I2004D . More recent 
spectroscopic monitoring campaigns for reverberation- 
mapping have been specifically designed to, among other 
things , recover velocity- de lay maps (e.g., iBentz et al.l 
I2010al iDennev et all 1201 ID . Thes e efforts have already 
been shown to be successful — IBentz et al.l ([2010bl ) 
show well-resolved velocity-delay map s for Arp 151 
using maximum-entropy methods, and iPancoast et al.l 
(!2012 f) used dynam ical modeling (first described by 
Panc oast et al.ll2011l ) to recover BLR structure informa- 
tion in Mrk 50. 

In late 2010, we carried out a four-month long rever- 
beration mapping campaign with the ultimate goal of 
recovering velocity- delay maps for several of the targets. 
Details of the data processing, light curves, mean time 
lags, and black hol e mass measuremen ts for these objects 
were published bv iGrier et al.l ([2012D . Here we present 
velocity-binned reverberation results for all five objects 
observed in our 2010 campaign and two-dimensional 
velocity-delay maps for four objects. All five targets 
and their coordinates and redshifts are listed in Table [TJ 
While our reverberation campaign was aimed primarily 
at investigating the H/3 emission line, in a few cases we 
recover velocity-delay maps for the H7 and He 11 A4686 
emission lines as well. 

2. DATA 

For this study we use spectroscopic and photometric 
data obtained during our reverberation campaign carried 
out at multiple institutions in late 2010. The primary set 
of spectra were obtained using the Boiler and Chivens 
CCD spectrograph on the MDM 1.3m McGraw-Hill tele- 
scope on Kitt Peak. We supplemented our spectroscopic 
continuum light curves with photometry from the 46- 
cm Centurion telescope at Wise Observatory of Tel- Aviv 
University and the 70-cm telescope at the Crimean As- 
trophysical Observatory (CrAO), and with spectra ob- 
tained from the 2.6m Shajn telescope at CrAO. All 
data were obtained between 2010 August 21 and 2011 
Ja nuary 7. T h e dat a processing is described in detail 
in IGrier et al.l (|2012f) . In short, we calibrated the re- 
duced spectra onto an absolute flux scale by assuming 
that the flux in the narrow [Oiii]A5007 emission line 
is constant. We created a reference spectrum for each 
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object using the spectra taken on photometric nights 
during the campaign. We then used the procedure of 
Ivan Groningen fc Wander j ()1992D to apply small wave- 
length shifts and scale each individual spectrum to match 
the flux in the [Oiii]A5007 emission line in the reference 
spectrum using a goodness-of fit estimator method. 
The resulting scaled spectra were used to create the driv- 
ing continuum light curves used throughout our analysis 
as well as the light curves used in our velocity-binned 
time series analysis. 

3. INITIAL H/3 VELOCITY-BINNED TIME SERIES 
ANALYSIS 

As a preliminary step, we first searched for gross kine- 
matic signatures by seeing if different parts of each emis- 
sion line show differen t time delays with r espect to the 
continuum. Following iDennev et al.l (|2009[ ). we divided 
the H/3 emission line into velocity bins as follows: First, 
we divide the H/3 emission lines in half at the zero- 
velocity line center (the systemic redshift), to separate 
red-shifted and blue-shifted signals. We then divide each 
line half into bins containing equal flux in the rms resid- 
ual spectrum, choosing the number of bins for each object 
based on the width of the rms line profiles. We create 
light curves for individual bins by integrating the flux 
within each bin in each scaled spectrum. The resulting 
light curves were analyz ed using JAVEL IN, the updated 
versio n of SPEAR (see iZu et all [20TTI and iGrier et al.l 
I2OI2I for details) , to measure the time delay in each light 
curve with respect to the continuum light curve. In 
essence, JAVELIN uses a statistical model of the con- 
tinuum light curve and its covariances and simultane- 
ously fits the continuum and line light curves assuming a 
simple top-hat transfer function. This leads to a statisti- 
cally well-defined means of interpolating the irregularly 
sampled line and continuum data that essentially aver- 
ages over all possible interpolated light curves weighted 
by their likelihood of fitting the data. In particular, the 
model fills gaps in the light curves in a well-defined man- 
ner with well-defined uncertainties. 

Our results from JAVELIN are shown in Figure [TJ The 
top panels show the rms residual line profile and the 
wavelength bins, and the bottom panel shows the dis- 
tribution of lags for each part of the emission line. As a 
check, we also ran the light curves for each bin through 
a cross correlation analysis routi ne tha t is widely used in 
reverberation mapping (e.g., Pe terson et al. 2004). The 
lags from the cross correlation were nearly identical to 
the ones obtained from JAVELIN. We show here only the 
JAVELIN results, since the JAVELIN uncertainties are 
much smaller than those obtained through cross correla- 
tion methods and the overall signatures are the same. In 
all five objects, we see velocity-dependent time lags in the 
H/3 emission line. Mrk 335, Mrk 1501, and PG 2130-1-099 
all show longer time lags on the blue side of the line pro- 
file than on the red side, which is a signature of inflowing 
gas. By contrast, 3C 120 shows a nearly symmetric pro- 
file, with small lags in the outer wings, and larger lags 
towards the line center, which is suggestive of a disk. 
Mrk 6 shows more complex structures that could be a 
co mbination of a disk plus i nfall, as previously suggested 
bv iDoroshenko et al.l ()2012[ ) for this object. 

4. TWO-DIMENSIONAL MEMECHO FITTING 



4.1. Preparation of Spectra 

To recover the velocity-delay maps, we used maximum 
entropy m ethods as impl emented in the pro gram MEME- 
CHO (see lHorne et al.lll99I and iHorndl 199 4 for details). 
To prepare the data for use with MEMECHO, we used 
software developed by Keith Horne called PREPSPEC. 
PREPSPEC applies corrections to account for wave- 
length shifts due to instrument flexure and differential 
refraction, spectrul blurring due to seeing and instrumen- 
tal resolution, and errors due to slit losses and changes 
in atmospheric transmission between epochs. To do this, 
we flt a calibration model to the spectra that accounts for 
both the spectral variability and the abovementioned sys- 
tematic errors. We model the spectra as the sum of flve 
different components: 1) a constant "mean" spectrum, 
which accounts for non-varying components of the con- 
tinuum and broad emission lines as well as starlight from 
the host galaxy, 2) non-varying narrow emission lines, 
3) a time-variable continuum, 4) time variable broad 
emission lines, and 5) factors to account for wavelength 
shift, spectral blurring, and scaling. After subtracting 
the models from the spectra, the resulting continuum- 
subtracted line proflle variations are used as inputs to 
MEMECHO. 

4.2. MEM Fitting 

In brief, MEMECHO flnds the "simplest" hnearized 
echo model that flts the observed continuum and 
emission-line spectral variations. This is accomplished 
by minimizing the function 

g = x'-2a^ (2) 

for a given model flt to a set of data. Here mea- 
sures the "badness-of-fit" between N data points with 
values Dk and the corresponding model predictions fik 
for those values, assuming Gaussian errors with known 
variances cr^ . The entropy S measures the "simplicity" of 
the model, elaborated below, and the regularization pa- 
rameter a controls the trade-off between these compet- 
ing requirements. MEMECHO adjusts a and the model 
parameters pi to achieve a user-specified value of ^ jN 
while maximizing S. 
In MEMECHO 's linearized echo model 

f'OO 

L{X,t) = Lo{X) + •^{X,T){C{t-r)-Co)dT, (3) 
Jo 

the parameters pi include the continuum light curve C(t) 
on an evenly spaced grid, the delay map at each wave- 
length \['(A,t), and the time- independent background 
spectrum io(A). Co is a reference continuum level, which 
we set at the median of the continuum light curve data. 
The entropy of the model is defined as 

S = ^Wi(j)i - Qi - piln{pi/qi)), (4) 

i 

where Wi are weights, pi are the positive parameters out- 
lined above, and qi are the default values of these parame- 
ters. Note that S is maximized when pi — qi. Minimizing 
Q gives 

= y^ "^'' 2'^'' + awj \n{pi/qi). (5) 
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Thus the model parameters pi are puUed by the data 
toward fik = and by the entropy toward pi = qi, with 
a adjusting the trade-off between the two. With defauh 
values Qi = {pi-iPi+iY^"^, the entropy "pulls" each pi 
toward the geometric mean of its neighbors, so that the 
entropy penalizes regions of high curvature and favors 
smooth functions with gaussian features and exponential 
tails. The weights Wi and default values qi are assigned 
with two parameters in MEMECHO: W to stiffen *(A, t) 
relative to C{t), and A to control the aspect ratio of 
features in ^(A, r). The A and W parameters are set by 
the user. Their qualitative effects on the velocity-delay 
maps are described below in Section fS. 2. II 

4.3. The Driving Continuum Light Curve 

One of the major practical issues with the use of 
MEMECHO is its tendency to introduce spurious fea- 
tures at gaps in the continuum light curve to drive ^ jN 
to the minimum value set by the user. This becomes a 
significant problem as the target /N decreases. Alter- 
ing MEMECHO parameters to stiffen (i.e., penalize rapid 
variations in) the driving light curve is sometimes help- 
ful, but we were still unable to produce velocity-delay 
maps with any discernible structure using the original 
continuum light curves. To provide stronger constraints 
on the driving light curve model, we ran our full con- 
tinuum light curves (including both spectroscopic and 
photometric data) through the JAVELIN modeling soft- 
ware of lZu et al.l ([2011). This method models the contin- 
uum with a damped random walk (DRW) model that has 
been demonstra ted to be a good statistical model of AGN 



variability (e.g. [Kelly et al.ll2009t iKoz lowski et "aLll2010l: 
MacLeod et al.l 

This allowed us to create highly sampled sim- 



ulated continuum light curves with reliable uncertain- 
ties that represent the range of the most likely contin- 
uum behavior within the gaps. Instead of providing 
MEMECHO with our original continuum light curve, 
we use our highly sampled JAVELIN mean light curve, 
which more strongly constrains the MEMECHO contin- 
uum model. This light curve is the probability-weighted 
mean of DRW light curves consistent with the data and 
the DRW model. Uncertainties on each point in the 
simulated light curve represent the standard deviation 
of probable light curves around this mean. These nar- 
row to match the measurement errors of the data points 
and then broaden as the gaps between the data become 
larger. The variability observed in the adopted mean 
simulated light curve is somewhat smoothed compared 
to an individual DRW model realization, but the error 
envelope applied to the mean light curve, representing 
the 1(7 deviations of the individual realizations about the 
mean, accounts for these differences. The original con- 
tinuum light curves and their corresponding simulated 
light curves are shown in Figure [5] 

Using the JAVELIN mean continuum light curve al- 
lows us to fit the variability of both the continuum and 
the rest of the spectrum to much higher degrees of ac- 
curacy and keep MEMECHO from introducing spurious 
features into our light curves and by extension, into the 
derived echo maps. Light curves, 1-dimensional MEME- 
CHO fits, and recovered delay maps at selected wave- 
lengths are shown in Figures |3][7] for all five objects. The 
wavelengths for the fits and delay maps shown in these 



figures were chosen to show the response in the red and 
blue wings as well as at the center of each emission line. 
One thing to note is that using the probability-weighted 
mean as the constraint on the MEMECHO construction 
of the continuum means that we are over-smoothing the 
light curves because we are neglecting the covariance of 
light curve deviations from the mean, which may weaken 
short-timescale variability. This will also cause some dif- 
ficulty in the fitting, as it can force the light curve mod- 
els to be smoother than the data. At the same time, the 
/N value reported by MEMECHO is no longer strictly 
valid — the deviations of any particular light curve model 
from the "mean" continuum light curve are highly cor- 
related, so, for example, having 10 consecutive points Icr 
from the light curve model is likely a la deviation, not a 
10(7 deviation. However, the velocity-delay maps do not 
change significantly if we use random individual DRW 
realizations of the light curve instead of the mean, indi- 
cating that the mean light curve can be used to produce 
accurate velocity- delay maps. While this solution is not 
ideal, it is currently the best method we have of dealing 
with gaps in the observed continuum light curves pend- 
ing a major effort to integrate the two distinct software 
packages (JAVELIN and MEMECHO). 

5. TWO-DIMENSIONAL VELOCITY-DELAY MAPS 

We were able to recover velocity-delay maps with 
MEMECHO for four out of the five galaxies observed. 
The best maps are for 3C 120 and Mrk 335, while the 
ones obtained for PG 2130-1-099 and Mrk 1501 are some- 
what less well-defined. The light curves for Mrk 6, on 
the other hand, do not seem to be well-fit by a simple 
echo model, and we were therefore unable to obtain two- 
dimensional velocity-delay maps for this object. Figure[6] 
shows the best model fits we were able to obtain for Mrk 
6. While there is some evidence for velocity-dependent 
structure in the one-dimensional delay maps shown, the 
models do not successfully fit any of the short-term vari- 
ations seen in this object. We were unable to improve the 
light curve fits by lowering the target x^/N, as MEME- 
CHO was unable to converge on a solution when we did 
so. 

The best-fit velocity-delay maps for the other four ob- 
jects, covering the full observed wavelength range for 
each, are shown in Figure [H The parameters used 
in MEMECHO to create the velocity-delay maps are 
given in Table [2J We also show more detailed velocity- 
delay maps for each individual emission line recovered for 
each object in Figures [91IT2] To aid the eye in compar- 
ing ionization-stratified BLR structure between different 
emission lines in each object, we also provide false-color 
maps in Figure 1131 with each color representing a differ- 
ent emission line as described in the captions. We also 
created simulated velocity-delay maps for a few different 
BLR kinematic models to compare qualitatively with our 
recovered velocity-delay maps. These simulated maps are 
shown in Figure 1141 ^-nd represent different BLR models 
around a black hole with Mbh = 1 x 10^ M0 . We see fea- 
tures reminiscent of these simple models in our recovered 
velocity-delay maps. 

5.1. Comments on Individual Objects 
5.1.1. Mrk 335 
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The velocity-binned analysis of Mrk 335 shows a defi- 
nite velocity-dependent lag signature (Figure [T]), and we 
clearly see similar structures in the velocity-delay maps 
of both the He ii and H/3 emission lines (Figure E]) . We 
see a chevron-shaped pattern in the Hen line, with a 
lack of prompt response in the center and shorter delays 
in the wings, which is consistent with the signatures ex- 
pected from an inclined disk or a spherical shell (Figure 
I14[) . In the H/3 emission, we see an asymmetric profile, 
with longer lags towards the blue end of the emission, and 
shorter lags towards the red end. This is suggestive of in- 
flowing g 0jS, 0jS demonstrated in our simple infall models 
given in Figure 1141 This also matches the signature we 
found in our initial velocity-binned analysis. The Heii 
emission is confined to smaller delays, as is expected from 
photoionization models of the BLR. Also consistent with 
disk structure is the response of He ii along a much wider 
velocity range than that of H/?. 

5.1.2. Mrk 1501 

Fewer kinematic details are apparent in the recovered 
maps for Mrk 1501 than for the other objects in our sam- 
ple. This may be a consequence of the noisier continuum 
and emission-line light curves obtained for this target. 
Nonetheless, we can still gain useful insights into the H/3 
and H7-emitting regions of the BLR in this object. As 
with Mrk 335, we see evidence in the velocity-delay maps 
for inflow, with longer lags in the blue and shorter lags to- 
wards the red end of the H/3 and H7 emission lines. This 
signature is closest to that of the "extended BLR" infall 
model shown in Figure [TH While this velocity-delay map 
is probably too low-resolution for any detailed modeling, 
we do see this same signature in our velocity-resolved 
analysis in Figure [TJ There is also evidence for radial 
stratification of the BLR — Figure [13] shows there is a 
much stronger H/3 response at longer time lags than is 
observed for the H7 emission. This is consistent with 
the idea that photoionization physics regulates the size of 
the BLR, as H7 is a higher-energy transition and would 
therefore need to be located closer to the ionizing source 
than H/3. 

5.1.3. 3C120 

Our velocity-delay maps for 3C 120, shown in full in 
Figure m are the cleanest maps we were able to recover. 
We see evidence for radial stratification, with He 11 show- 
ing the shortest delays and H7 and H/3 emitted at pro- 
gressively larger radii. Figure [T3] shows all three emis- 
sion lines on one scale, with H/3, H7, and He li shown in 
red, blue, and green, respectively, highlighting the radial 
stratification. The lack of a prompt response at line cen- 
ter in both Balmer lines indicates a deficit of material 
along the line of sight. The shape of the hydrogen re- 
sponse is consistent with signatures expected from an in- 
clined disk or a spherical shell (Figure [Til). This is similar 
to th e signal found in NG C5548 (Hornc ct al. 1991 ), Arp 
151 ([Bentz et all l2010b( l. and Mrk 50 (Pancoa st et al.l 
I20TI . We also see an asymmetry in the strength of the 
response in both the He 11 and H/3 emission profiles, nom- 
inally indicating inflow. The smaller velocity range of the 
H/3 and H7 emission with respect to Hell, combined with 
the radial stratiflcation signatures, is consistent with disk 
structure. 



5.1.4. PG 2130+099 

The velocity-delay map for PG 2130-1-099 is somewhat 
noisier than that of 3C 120, and the model fits are not 
as good. However, the structure of the map is worth 
noting. While reverberation delays for the He 11 emission 
remain unresolved for this data, both the H/3 and H7 
emission lines clearly show velocity-resolved delay struc- 
ture, with strong asymmetries as a function of velocity 
across the emission-line profile. This is similar to the 
H/3 emission we see in Mrk 335 and the Hen emission 
in 3C 120. This asymmetry, with longer lags at the blue 
end of the emission line and shorter lags to the red, is 
suggestive of inflowing gas, and matches the "Infall (less- 
extended BLR)" model in Figure [14] quite well. In our 
analysis of the velocity-binned sections of the H/3 emis- 
sion line (Figure [T]), we see this same structure, so while 
the delay maps for this object are probably not good 
enough for any detailed modeling, they are consistent 
with the signatures of inflowing material already seen in 
this object. 

PG 2130-1-099 h as long been a curiosity. Early on, 
iKaspi et al.l (|2000f ) measured a time lag of ^180 days in 
this object, placing it well above the R bl^- L relation- 
ship. Later studies bv lGrier et al.l (|2008D and lGrier et al.l 
(|2012|) found much shor ter lags on the o rder of tens 
of days: Most recently, iGrier et all (|2012[ ) reported a 
mean H/3 time delay of 12.8to'9 days. These stud- 
ies attributed the discrepancies to undersampled light 
curves combined with long-term sec ular changes in th e 
H/3 equivalent width in the data from lKaspi et al.l (|2000f) . 
However, even with the new, shorter lag measurements, 
PG 2130-1-099 is still a major outlier from the i?BLR-i 
relation, as it is n ow positioned far below the relation 
(iGrier et al.l l2012t ) . Despite the higher sampling rate 
of the more recent campaigns, ambiguities remain as to 
whether the measured H/3 lags represent the true mean 
BLR radius, as the light curves were missing data at key 
points in time. We see in our velocity-delay map that 
the majority of the response in the H/3 emission seems to 
be centered on a delay of ~30 days (Figures [12] and [13]) . 

To investigate this, we ran a one-dimensional delay 
map analysis of PG2130-h099 in MEMECHO to look 
for an indication of where the true lag lies. Figure [15] 
shows the model continuum light curve env elope in the 
bottom panel, and the H/3 light curve from IGrier et al.l 
(2012) in the top right panel; the top left panel shows the 
delay map recovered by MEMECHO. The MEMECHO 
model fits the data fairly well in this case, and there are 
two clear peaks in the delay map. The stronger peak is 
centered around 12.5 days, and the slightly weaker peak 
is centered at 31 days. We compare this with the two- 
dimensional velocity-delay map (Figures [H] [121 and [13]) , 
which shows a large signal on the blue side of the emission 
concentrated at around 30 days and a fainter signal to the 
redward side that stretches down to shorter lags. A plau- 
sible model reproducing these results is a nearly face-on 
disk with the emitting gas located at around 30 light- 
days, combined with a strong infiowing gas component 
not necessarily within the plane of the disk. Including 
an inflow signature when we measure the flux of the en- 
tire H/3 emission line could cause our result to be skewed 
towards shorter mean lags, when the true distance of the 
virialized gas is closer to ~30 days. Because of the lower 
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quality and coarser sampling of the light curves for this 
object, we will likely be unable to model this structure 
in much more detail. However, it is clear from the delay 
map that the majority of the H/3 signal comes from a 
radius of --31 hght-days. This radius puts PG2130+099 
much closer to the i?BLR~i relation. This also increases 
the black hole mass estimate for PC 2130+099 by a fac- 
tor of about 2.4, putting it at about 10^ Mq. This would 
place PG 2130+099 definitely within the expected scatter 
of the Mbh-ct* relation. 

5.2. MEMECHO Parameters and Settings 
5.2.1. The A and W Parameters 

MEMECHO gives the user control over various aspects 
of the fitting process. We first consider the weights Wi 
(see EquationHl), which are implemented in MEMECHO 
as the user-controlled W and A parameters. The W pa- 
rameter controls the weight given to pixels in ^(A, t) rel- 
ative to the weight given to the pixels in the continuum 
model C{t). Increasing W makes the delay map stiffer 
and allows for more flexibility in the continuum model. 
A affects the aspect ratio of features in the velocity-delay 
map \E'(A, r). The default values for ^'(A,t) are geomet- 
ric means of neighboring pixels, with A increasing the 
weight of neighbor pixels in the A direction relative to 
those in the t direction. 

It is important that the structures we see in our de- 
lay maps are not dependent on our choices of A and 
W . To verify that our results are robust, we produced 
MEMECHO models altering the A and W parameters 
and changing the wavelength and velocity grids of the 
transfer function models. The panels of Figure [16] show 
the resulting velocity-delay maps for 3C 120 as we vary 
these parameters. For the most part, the general shape 
of our velocity-delay map is not affected by changes in 
A or W . However, if A is made too large, we start to 
over-fit the data and introduce spurious features into the 
maps. The stability with respect to the W parameter is 
expected, since the continuum light curve we are using 
is very highly sampled, and therefore very highly con- 
strained. Since these behaviors are generic to all four 
objects, we only show the results for 3C 120. 

5.2.2. The Mimmum x^/N 

We also investigated the effect of varying the target 
X^/N of the MEMECHO solution, as the final resolution 
of the recovered delay maps is controlled by this parame- 
ter. The choice of the target /N is a trade off between 
smoother, lower delay resolution maps for higher x^ jN , 
and sharper but less reliable structure in the maps at 
lower x^ /N . Because we are using over-sampled contin- 
uum model light curves from JAVELIN, the actual value 
of x^ /N is no longer strictly valid, because the continuum 
light curve data points are not fully independent. There- 
fore, our chosen values oi x'^/N reflect the best trade-off 
between delay-map resolution and reliable structure in 
the fits that we could obtain. Figure [17] shows velocity- 
delay maps for 3C 120 as we vary the x^ /N , from left to 
right. As we lower the target x^ /N , the structure even- 
tually becomes more complex but also less reliable, as 
the fits are now producing structures to model the noise 
in the line light curves. However, the basic structure of 
the velocity-delay maps is robust to reasonable changes 



in the target x^ /N . 

5.2.3. Minimum Allowed Time Delays 

When we vary the minimum lag Tmin that we allow 
MEMECHO to consider, we also see a trend worth not- 
ing. Figure [IT] shows that when Tmin is set to zero (top 
panels), the velocity- delay structure in the map also ex- 
tends all the way to zero. However, when negative lags 
are allowed (middle and bottom panels), we see that the 
structure in the velocity-delay maps, particularly when 
considering the chevron-shaped structure in the H/3 and 
H7 emission lines, peaks at a positive delay and is lower 
at T = but does not extend to 0. This has to do with 
the way MEMECHO deals with its delay map models. 
The default is set to "pull" the response down to zero 
at the ends of the delay map when Tmin and r^ax are 
not equal to zero. Therefore, when Tmin — 0, the default 
cannot pull down at the low-r end, so the lowest-entropy 
map is an exponential function of r. With Tmin < 0, the 
entropy pulls down on both ends of the delay map, which 
then creates a Gaussian peak at positive delays. 

Because of this, different values of r^-m can affect the fi- 
nal velocity-delay maps produced, as this effectively reg- 
ulates the behavior of the model on some level when try- 
ing to smooth the delay maps to the default image. We 
currently have no means to evaluate which value of Tmin 
results in the "correct" velocity-delay map — it is there- 
fore worth considering velocity-delay maps with various 
■''min parameters, as we do here, when drawing conclu- 
sions about the BLR structure signatures, and any de- 
tailed modeling and interpretations should consider and 
evaluate the differences in these maps. In this work, we 
make only qualitative statements regarding the signa- 
tures seen in the velocity-delay maps. Whether or not 
the response actually reaches zero at any point in the 
3C 120 maps is unknown, but for our purposes it is reas- 
suring that all of the maps show the same basic structures 
and kinematic signatures. From Figure 1171 comparing 
results for Tnii„ — —10,-1, and 0, we can also evaluate 
the evidence for a deficit of prompt response from the 
resulting delay maps created with each value of Tmin . At 
all three values of Tmin, we see Gaussian delay distribu- 
tions at positive lags in the core of H/3 and H7 in 3C 120, 
but not for He li — thus the data provide evidence for a 
deficit of prompt response in the center of H/3 and H7, 
but not for He 11. This behavior is consistent with our 
expectations for the stratification of the BLR. 

6. SUMMARY AND CONCLUSIONS 

We have presented velocity-delay maps constructed 
from the line and continuum variations observed in four 
objects from our 2010 reverberation campaign. These 
maps provide new insights into the structure of the BLR 
and constitute a dramatic increase in the number of ob- 
jects that have at least some information on the velocity 
field of t he BLR. Al ong with the velocity-delay maps for 
Arp 151 (iBentz et a l. 2010b) and the models of Mrk 50 
(jPancoast et al.ll20r2) , these velocity- delay maps provide 
the strongest constraints on the structure of the BLR. 
Our velocity-delay map for 3C 120 shows very similar 
structure to the map of Arp 151, which also shows a lack 
of prompt response in the Balmer lines. The continuum 
light cu rves for Arp 151 were well-enough sampled that 
IBentz et al. (,2010bl were able to recover a velocity-delay 
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map using the original continuum light curve instead of 
substituting in simulated ones. The similarities between 
the maps for Arp 151 and our maps lends weight to the 
reliability of our results using the simulated continuum 
light curves. We also see very asymmetric profiles in both 
the Balmer and high-ionization emission that is sugges- 
tive of infalling gas in all of our objects. In 3C 120 and 
Mrk 335, the transfer function structure for the Balmer 
lines differs from that for the HeilA4686 emission line, 
suggesting different structures dominating at different 
BLR radii. In all cases where our data are of sufficient 
quality to constrain the structure of the BLR, we see 
clear evidence of infall and rotation, both of which result 
from the gravitational influence of the black hole. As 
gravitationally dominated motion is the key assumption 
of reverberation mapping, our results strongly support 
the reliability of black hole mass estimates derived from 
reverberation mapping. Detailed modeling of our most 
well-defined velocity- delay maps and a complete discus- 
sion of the implications will follow in a future work. 
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TABLE 1 

Object List 



Object 




RA 


DEC 


z 




{J2000) 


(J2000) 




Mrk 335 


00 


06 19.5 


+20 12 10 


0.0258 


Mrk 1501 


00 


10 31.0 


+ 10 58 30 


0.0893 


3C 120 


04 


33 11.1 


+05 21 16 


0.0330 


Mrk 6 


06 


52 12.2 


+74 25 37 


0.0188 


PG 2130+099 


21 


32 27.8 


+10 08 19 


0.0630 



TABLE 2 
MEMECHO Parameters'' 



Object 


A 


W 






Mrk 335 


0.1 


1.0 


-1 


2.8 


Mrk 1501 


1.0 


1.0 





1.8 


3C120 


1.0 


1.0 





2.4 


PG 2130+099 


0.1 


1.0 





2.1 



As defined in Section |4] 
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Fig. 1. — Velocity-binned reverberation lag results. The top panels show the rms residual spectrum for each object, with the edges of the 
bins designated by vertical dashed lines. The bottom panels show the mean time delays measured for each bin. T he zero-velocity ce nter of 
the H/9 emission line is shown by the dotted line. The horizontal solid line shows the average time lag reported in lGrier et al.l gm2), with 
uncertainties as horizontal dashed lines. Error bars in the velocity direction show the width of the bins. 
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HJD (-2450000) 

Fig. 2. — Original continuum light curves (gray data points) from [Grier et all pOl^ *) and the JAVELIN mean continuum model (black 
solid line) used in the MEMECHO analysis. The dotted line shows the standard deviation about the mean light curve from individual 
model realizations. Fluxes are given in units of 10~^^ erg s~^ cm~^A~^. 
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Fig. 3. — Best MEMECHO fits to the spectra of Mrk 335. The left panels show the delay maps at each selected rcst-framc wavelength, 
given in the top righthand corner of each panel. The right panels show the emission-line light curve at these selected wavelengths and the 
MEMECHO fit to the light curve. The X of tiic fit for each light curve is shown in each panel. The bottom panel shows the original 
continuum light curve (black error bars), the error envelope from the simulated light curve showing the standard deviation about the mean 
(gray envelope), and the continuum model from MEMECHO (solid black line). 
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Fig. 4. — Best MEMECHO fits to the spectra of Mrk 1501. Panels and symbols are the same as in Figure[3] 
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Fig. 5. — Best MEMECHO fits to the spectra of 3C 120. Panels and symbols are the same as in Figure[3] 
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Fig. 6. — Best MEMECHO fits to the spectra of Mrk 6. Panels and symbols are the same as in Figure [3] 
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Fig. 7.— Best MEMECHO fits to the spectra of PG 2130+099. Panels and symbols are the same as in Figure|3l 
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Fig. 8. — Our best-fit velocity-delay maps over the full rest-frame wavelength range for each object (grayscale), projections onto the 
rest-frame wavelength axis (bottom panel) and the time-delay axis (right panel). 'I' (A) is the overall response added up at each wavelength, 
and ^(r) is the overall response of all emission lines added together at a given r. 
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Fig. 9. — Velocity-delay maps for bot h emission lines se en in Mrk 335. The dotted lines show the "virial envelope", V^tc/G = 4.6 X 10^ 
M0, measured from the mean time lag l IGrier et aI.II20l"3 ). 
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Fig. 10. — Velocity-delay maps for both e mission lines se en in Mrk 1501. The dotted lines show the "virial envelope" , V'^tc/G = 3.3 X 10^ 
M0, measured from the mean time lag ((Grier et al.ll20l"3 '). 
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Fig. 11. — Velocity-delay maps for all three emission lines seen in 3C 120. The dotted lines show the "virial envelope", V^tc/G = 1.2 X 10^ 
Mq, measured from the mean time lag (Grier et al. 2012). The velocity-delay map for H7 is truncated at V'~3250 km/s because the 
wavelength range of our spectrograph did not cover this far into the red. 
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Fig. 12. — Velocity-delay maps for all three emission lines seen in PG 2130+099. The dotted lines show the "virial envelope", V^rc/G = 
8.3 X 10^ Mq , measured from the mean time lag IIGrier et al.|[20T3) . 
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Fig. 13. — False-color velocity-tielay maps for Mrk 335, Mrk 1501, 3C 120, anti PG 2130-1-099. The dotted lines in each panel correspond 
to virial envelopes for each object as listed in Figures l9ll2l The H/3 emission is shown in red, H7 emission in green, and Hell A4686 emission 
in blue. 
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Fig. 14. — Velocity-delay maps for simple BLR models of H/3 emission around a 1x10^ Mq black hole. Panel (a) shows highly beamed 
emission (typical for H/3, see Ferland et al. 1992) from gas in free-fall motion. The infalling gas is distributed in a spherical shell with inner 
and outer radii of 15 and 20 light-days, inclined at an angle of 45 degrees. Panel (b) shows the same infall model but with a more-extended 
BLR with inner and outer radii of 5 and 20 light-days, respectively. Panel (c) shows a map for an edge-on Keplerian disk with inner and 
outer radii of 5 and 20 light-days, and panel (d) shows a map for a fully illuminated thick spherical shell of Keplerian circular orbits, with 
inner and outer BLR radii of 5 and 20 light-days. 
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Fig. 15. — One-dimensional delay map for PG 2130-(-099. The bottom panel shows the simulated continuum light curve used in the 
MEMEC HO analysis , with the errors shown as the black envelope. The top right panel shows the light curve for the entire H/3 emission 
line from lGrier et ahl II2012I) . The top left panel shows the one-dimensional delay map from MEMECHO. 
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Fig. 16. — Velocity-delay maps for 3C 120. Each panel shows the best map found as we vary A and W by factors of 100. Increasing A 
smooths the maps more strongly in A and less in r, while increasing W reduces the importance of fitting the continuum model relative to 
the delay maps. All models are converged to the same overall goodness-of-fit fN . Our adopted velocity-delay map is that of the middle 
panel. 




Fig. 17. — Velocity-delay maps for 3C 120, varying the degree of smoothing and the minimum lag allowed in the model (rmin)- Increasing 
values of /N mean that the entropy term is more heavily smoothing the map. Our adopted velocity-delay map is that of the top middle 
panel. 



